Quantum back-action in measurements of zero-point mechanical oscillations 
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Measurement-induced back action, a direct consequence of the Heisenberg Uncertainty Principle, is the defin- 
ing feature of quantum measurements. We use quantum measurement theory to analyze the recent experiment 
of Safavi-Naeini et al. [Phys. Rev. Lett. 108, 033602 (2012)], and show that results of this experiment not 
only characterize the zero-point fluctuation of a near-ground-state nanomechanical oscillator, but also demon- 
strate the existence of quantum back-action noise — through correlations that exist between sensing noise and 
back-action noise. These correlations arise from the quantum coherence between the mechanical oscillator and 
the measuring device, which build up during the measurement process, and are key to improving sensitivities 
beyond the Standard Quantum Limit. 

PACS numbers: 



I. INTRODUCTION 

Quantum mechanics dictates that no matter or field can stay 
absolutely at rest, even at the ground state, for which energy 
is at minimum. A starting point for deducing this inevitable 
fluctuation is to write down the Heisenberg Uncertainty Prin- 
ciple 



which leads to: 



[x,p] =ih, 



Ax-Ap > ft/2. 



(1) 



(2) 



Here x and p are the position and momentum operators, while 
Ax and Ap are standard deviations of position and momen- 
tum for an arbitrary quantum state. Eq. |2]) means we cannot 
specify the position and momentum of a harmonic oscillator 
simultaneously, as a point in classical phase space — the os- 
cillator must at least occupy ft/2 area in the phase space. If 
the oscillator has mass of m and eigenfrequency of co m , then 
in the Heisenberg picture we can write 



[x q (t),x q {t')] = 



which leads to: 



Ax q (t) ■ Ax q (t') > 



ih sin a> m (t' — t) 
mco m 



ft|sina)„,(f'-f)| 

ImCOm 



(3) 



(4) 



with x q (t) being the Heisenberg operator of the oscillator po- 
sition, quantum-mechanically evolving under the free Hamil- 
tonian. Here Ax q (t) is the standard deviation of x q (t) for an 
arbitrary quantum state. Eq. Q means the position of a freely 
evolving quantum harmonic oscillator cannot continuously as- 
sume precise values, but instead, must fluctuate. This fluctua- 
tion carries the zero-point mechanical energy of hco m /2. 

As a key feature of quantum mechanics, zero-point fluctua- 
tion of displacement is an important effect to verify when we 
bring macroscopic mechanical degrees of freedom into their 



ground states HKH1- Needless to say, a continuous observa- 
tion of the zero-point fluctuation of a macroscopic mechanical 
oscillator requires superb displacement sensitivity. 

However, what constitutes an "observation of the quantum 
zero-point fluctuation" is conceptually subtle. Eqs. <|3j and 
Q, which argue for the inevitability of the zero-point fluc- 
tuation, also dictate that the "exact amount" of the zero-point 
fluctuation cannot be determined precisely. More specifically, 
if we use a linear measurement device to probe the zero-point 
fluctuation, which has an output field of y(t ), then we must at 
least have 



[y(t)J(t')}=0 



(5) 



at all times, in order for y(t) to be able to represent experi- 
mental data string — with measurement noise simply due to 
the projection of the device's quantum state into simultaneous 
eigenstates of all {y(t) : t € R}. This means y must be written 
as 



y{t) = e{t) - 



-x q {t) 



(6) 



with non-vanishing additional noise (error) e(f), which con- 
sists of degrees of freedom of the measurement device and 
compensates the non-vanishing commutator of x q In addi- 
tion, during the measurement process, actual evolution of the 
mechanical displacement x must differ from its free evolution 



x q . This is because 



ih 



(7) 



is also the classical response function of x to an external force: 
any device that attempts to measure x by coupling it with an 



* We note that Ozawa has developed a different formalism to quantify the is- 
sues that arise when attempts are made to measure non-commuting observ- 
ables like x q (t) (9l ll01 . However, we have chosen to adopt the Braginsky- 
Khalili approach [lj], because it is immediately applicable when the non- 
commuting observable is acting as a probe for an external classical force. 
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external observable F, which introduces a term proportional 
to xF into the Hamiltonian, will have to cause non-zero dis- 
turbance. For this reason, we can expand the measurement er- 
ror e into two parts: £ — the sensing noise that is independent 
from mechanical motion and xba — additional disturbance to 
the mechanical motion from the measurement-induced back- 
action, and rewrite y(t) as: 



y(t) = z{t) +*ba(*) +x q {t) = z(t) + x{t) . 



(8) 



The mechanical displacement under measurement is therefore 
a sum of the freely-evolving operator x q plus the disturbance 
xba due to back action noise, namely, x(t) = x q (t) + jcba(0- 

The above lines of reasoning lie very much at the heart 
of linear quantum measurement theory, pioneered by Bra- 
ginsky in the late 1960s aiming at describing resonant-bar 
gravitational -wave detectors ifTTl fl2l . and later adapted to the 
analysis of laser interferometer gravitational-wave detectors 
by Caves lfl3l . A key concept in linear quantum measure- 
ment theory is the trade-off between sensing noise and back- 
action noise, which gives rise to the so-called Standard Quan- 
tum Limit (SQL). For optomechanical devices, sensing noise 
takes the form of quantum shot noise due to discreteness of 
photons, while the quantum back-action is enforced by quan- 
tum fluctuations in the radiation pressure acting on the me- 
chanical oscillators |[T3l , which is therefore also called quan- 
tum radiation-pressure noise. It has been shown that the SQL, 
although not a strict limit for sensitivity, can only be surpassed 
by carefully designed linear measurement devices which take 
advantage of quantum correlations between the sensing noise 
and the back-action noise. 

Observing signatures of quantum back-action, achieving 
and surpassing the associated SQL in mechanical systems are 
of great importance for the future of quantum-limited metrol- 
ogy, e.g., gravitational-wave detections IIT41 - I221 . At the mo- 
ment, it is still experimentally challenging to directly observe 
quantum radiation-pressure noise in optomechanical devices 
due to high levels of environmental thermal fluctuations, and 
there are significant efforts being made toward this RH8l l23l . 
One approach proposed by Verlot et al. (U is, instead, to 
probe the quantum correlation between the shot noise and 
the radiation-pressure noise, which, in principle, is totally im- 
mune to thermal fluctuations. 

In this article, we analyze a recent experiment performed 
by Safavi-Naeini et al. 11241 . in which a radiation-pressure- 
cooled nanomechanical oscillator — the movable mirror of 
a high-finesse cavity — is probed by a second beam of light, 
detuned from the cavity, for its zero-point mechanical oscilla- 
tion. The output power spectrum of the second beam, near the 
mechanical resonant frequency, serves as an indicator of the 
oscillator's zero-point motion. It was experimentally observed 
that when the second beam is detuned on opposite sides from 
the cavity resonance, the output power spectra turn out to be 
different. Using theory of linear quantum measurements, we 
will show that this experiment not only probes the zero-point 
fluctuation of the mechanical oscillator at nearly ground state, 
but also illustrates vividly the non-trivial correlations between 
sensing noise and back-action noise — an much sought-after 



effect in the gravitational-wave-detection community. Its con- 
tribution to tlie output spectrum is equal to the zero-point fluc- 
tuation for one detuning of the readout beam, and exactly op- 
posite for the other detuning. 

The outline of this article goes as follows: in Sec.|ll] we will 
give a brief overview of the experiment by Safavi-Naeini et 
al, and present an analysis of this experiment using quantum 



measurement theory; in Sec. Ill we will more broadly discuss 
the nature of mechanical zero-point fluctuation, show that in 
attempts to measure the zero-point fluctuation, the contribu- 
tions from sensing-back-action noise correlations can generi- 
cally be comparable to the zero-point fluctuation itself. In ad- 
dition, we will discuss linear quantum measurement devices 
that use a near-ground-state mechanical oscillator as a probe 
for external classical forces near its resonant frequency, and 
show the limitation on the measurement sensitivity imposed 
by the zero-point fluctuation and the connection to the SQL; 
we will conclude in Sec.lIVI 



II. 



A TWO-BEAM EXPERIMENT THAT MEASURES 
ZERO-POINT MECHANICAL OSCILLATION 



In this section, we describe in Sec. |II A| the experiment 
performed by Safavi-Naeini et al., put its results into the 



framework of linear quantum measurement theory in Sec . II B 



and provide a detailed analysis in Sec. II C In Sec. II D 



we 



will comment on the connection between the viewpoint from 
quantum measurement and the scattering picture presented in 
Ref.l24l. 



A. Experimental setup and results 

In the experiment, two spatial optical modes are cou- 
pled to a mechanical vibrational mode in a patterned sili- 
con nanobeam. One spatial mode — the cooling mode — 
is pumped with a relatively high power at a "red" detun- 
ing (lower than resonance), and is used to cool the mechan- 
ical mode via radiation pressure damping |2|; the other cav- 
ity mode — the readout mode — has a much lower power 
and is used for probing the mechanical motion. The readout 
laser frequency ©/, is detuned from the resonant frequency (£>,- 
of the readout mode by either +(Q,„ or — ©„,. The observed 
spectra of the readout laser are asymmetric with respect to the 




FIG. 1: Figure illustrating the observed spectra of the readout laser 
in the positive-detuning case (left) and the negative-detuning case 
(right). 
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detuning A = co r — ©/,-. Specifically, in the positive-detuning 
case — A = co m , the spectrum has a smaller amplitude than 
that in the negative-detuning case. The area I + enclosed by 
the spectrum in the positive-detuning case, after subtracting 
out the noise floor away from the mechanical resonant fre- 
quency, is proportional to the thermal occupation number («) 
of the mechanical oscillator, while, in the negative-detuning 
case, the enclosed area is /_ °c (n) + 1. Such asymmetry is 
illustrated in Fig.[T] In Ref. Il24l . we introduced the following 
figure of merit to quantify the asymmetry: 



1 = T^ 1 = T~\ 
I, (n) 



(9) 



We interpreted this asymmetry as arising from the quan- 
tized motion of the mechanical oscillator. The asymmetry is 
thus assigned to the difference between the phonon absorption 
rate, proportional (n), and the emission rate, proportional to 
(n) + 1 . This is completely analogous to that used for calibra- 
tion of motional thermometry of ions/atoms trapped in electri- 
cal/optical traps 125H281 . Additionally, these scattering pro- 
cesses have an underlying physics similar to bulk nonlinear 
Raman scattering processes used in spectroscopic analysis of 
crystals [29, 30 1, where an ensemble of vibrational degrees of 
freedom internal to the molecular structure of the system in- 
teracts with incident light. Typically in these nonlinear optics 
experiments, photon counters are used to keep track of the 
(anti-)Stokes photons. In contrast, in our experiment, a het- 
erodyne measurement scheme was used, to find the amplitude 
quadrature of the readout mode. Interestingly, by choosing the 
detuning A = ±0) m and in the resolved-sideband regime, spec- 
tra of the amplitude quadrature are equal to emission spectra 
of the (anti-)Stokes photons plus a constant noise floor due 
to vacuum fluctuation of the light — the shot noise. We will 



elaborate on this point in Sec. II D and show explicitly such 



a connection. Intuitively, one can view the cavity mode as an 
optical filter to selectively measure the emission spectra — for 
A = a),,,, the anti-Stokes process is significantly enhanced as 
the emitted photon is on resonance with respect to the cavity 
mode, and one therefore measures the spectrum for the anti- 
Stoke photons; while for A = — ©„,, the spectrum of the Stokes 
photon is measured. 



B. Interpretation in terms of quantum measurement 



Here we provide an alternative viewpoint to Ref. [24|, em- 
phasizing on the role of quantum back-action and its relation 
to quantization of the mechanical oscillator. First of all, we 
separate the experimental system into two parts. The cooling 
mode, the mechanical oscillator, and the environmental ther- 
mal bath the oscillator couples to (the left and middle boxes in 
Fig. |2]i together is the first part, which can be viewed as pro- 
viding an effective mechanical oscillator nearly at the ground 
state, but with a quality factor significantly lower than the in- 
trinsic quality factor of the mechanical mode. It is the zero- 
point fluctuation of this effective oscillator that we shall be 
probing. The second part of the system consists of the readout 
mode (the box on the right of Fig. EJ, which couples to the 
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FIG. 2: Figure illustrating the relation among different parts of the 
optomechanical system in the experiment. The thermal heat bath and 
the cooling mode together create an effective quantum heat bath for 
the mechanical oscillator which in turn couples to the readout mode. 



effective oscillator (the first part of the system) through dis- 
placement i alone. The second part provides us with an output 
y, which contains information about the zero-point fluctuation 
of the effective mechanical oscillator. 



1. The Mechanical Oscillator Near Ground State 

Let us focus on the first part of the system (left two boxes 
of Fig. [2l, the effective mechanical oscillator (Because this 
will be a stand-alone subject of study in later discussions, we 
shall often ignore the word "effective"). The environmental 
heat bath and the cooling mode together form a quantum heat 
bath with fluctuation close to the zero-point value. In steady 
state, the "free" mechanical displacement is determined by its 
coupling to this bath ("free" means in absence of the readout 
mode): 



X{t-t')F q {t')At' . 



(10) 



Here % is the response function of the mechanical oscillator: 



mxt')] 



-K m \t-t'\/2 



ih 



sino) m (f — t 1 ) 
mco m 



(11) 



Note that we have an additional decay factor compared with 
Eq. ([3]) which describes an idealized free oscillator. The de- 
cay rate K m here is determined jointly by the intrinsic decay 
rate of the mechanical mode, and the optomechanical interac- 
tion between the mechanical mode and the cooling mode. The 
force F q lumps together the fluctuating force acting on the me- 
chanical mode by the environmental heat bath and the cooling 
mode. If the oscillator approaches the ground state only af- 
ter applying the cooling mode, then one can show that F q is 
dominated by fluctuation of the cooling mode. 

The above two equations state that for a realistic mechani- 
cal oscillator with non-zero decay rate, its zero-point fluctua- 
tion in the steady state can be viewed as driven by the quantum 
heat bath surrounding it. We will returning to this prominent 



feature of linear quantum systems later in Sec. Ill A 



2. The Quantum-Measurement Process 

Let us now move on to the second part of the system (right 
box of Fig.|2]i, in which the readout mode serves as a linear po- 
sition meter that measures the mechanical displacement. We 
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can rewrite the disturbance xba in Eq. ([H} in terms of the back- 
action force Fba arising from radiation-pressure fluctuation of 
the readout mode, namely, 



*ba(0 



X(t-t')F BA (t')dt'. 



(12) 



We have assumed that the readout mode does not modify the 
dynamics of the oscillator, which is a good approximation for 
the low pumping power used in the experiment. Written in the 
frequency domain, the readout mode output y [cf. Eq.([8]l] is 



where 



y(co) ^ z(co) + x((o)Fba{co) + x{(o)F q {co). 



X(a>) 



m(co 2 - 0)2 + iK m C0) 



(14) 



is the Fourier transform of &(t)x(t), with © the Heaviside 
function, i.e., the positive half of x(t ) (even though ^(f) exists 
for both t > and t < 0). The spectral density S yy (co) of y then 
reads: 



Syy = S zz + 2Re[ X *S zF ] + \ X \ 2 S B FF + |*| 2 ^ 



(15) 



Here these single-sided spectral densities are defined in a sym- 
metrized way (see Appendix [A]), which guarantees bilinearity 
for the cross spectrum and positivity for self spectrum; S zz and 
Spp are the sensing-noise and back-action force noise spec- 
trum, respectively; S z f is the cross correlation between z and 
/ba; the force spectrum of the effective quantum heat bath 
made up by the environmental heat bath and the cooling mode 
is given by: 



(4(n) +2)hmK m C0„, 



(16) 



and (n) is the thermal occupation number. 



3. Asymmetry between Spectra 

Experimentally, it was observed that the output spectra Syy 
for the two opposite detunings, A = ±0),,,, are different — 
given the same thermal occupation number for the oscillator, 



M®)lA=-<»,„ ¥*Syy((0)\ 



A=Wm 



(17) 



As we will show in the Sec. |IIC| that follows, when we flip the 
sign of the detuning A of the readout beam, the only term in 
Syy that changes is S z f — the correlation between the sensing 
noise and the back-action noise. According to Eq. ( f34| >, we 
have 



S zF ((0): 



-ih 



a> 

A' 



(18) 



in the resolved-sideband regime with the cavity bandwidth 
K r <C CO,,,, which is the case in the experiment. The asymmetry 
factor defined in Eq.|9]) is given by: 



2/Re[x*(^-S+)]d£0 



J [\x\ 2 S c > F +2Re(x*St F )]dco 



(19) 



Here Si. is defined by Sf F = S z f |a=±<h„> an d i n particular 
around the mechanical resonant frequency *»„,, where Szp con- 
tribute to the above integral, 



(20) 



The asymmetry, or effect of quantum correlation S z f, is 
most prominent when the thermal occupation number ap- 
proaches zero. Indeed, if we focus on the quantum fluctuation 
by taking the limit of (n) — > 0, we obtain 



(13) j 2Re( X * Sf F )dco 



\x\ z s 



FF\{n)=0 



=n dO). 



(21) 



In other words, at the quantum ground state, contribution of 
the quantum correlation S z f to the readout spectrum Syy has 
the same magnitude as that of the zero-point fluctuation, while 
the sign of the correlation term depends on the sign of the de- 
tuning of the readout beam. This means not only has the ex- 
periment probed the zero-point fluctuation of the mechanical 
oscillator, it has also demonstrated non-trivial correlations 
between sensing noise and back-action noise at the quantum 
level. 



C. Detailed theoretical analysis 

In this section, we supply a detailed calculation of the quan- 
tum dynamics and the output spectrum of the experiment. The 
dynamics for a typical linear optomechanical device has been 
studied extensively in the literature ll3T1433l : however, they 
have been focusing on quantum state of the mechanical os- 
cillator in ground-state cooling experiments, instead of treat- 
ing the optomechanical device as a measurement device. Here 
we will follow Ref. |34] and derive the corresponding input- 
output relation — the analysis is the same as the one of quan- 
tum noise in a detuned signal-recycling laser interferometer 
which can be mapped into a detuned cavity fl31 [35l l36l . We 
will focus only on the interaction between the readout cavity 
mode and the mechanical oscillator — the cooling mode and 
the thermal heat bath is taken into account by the effective 
dynamics of the oscillator as mentioned earlier. 

The Hamiltonian of our optomechanical system can be 
written as IT3TH331 : 



M 3 — h co r d T d + J^ Kr + h Goxa'a 



El. 

2m 



2 «2 

mfflj - 



(22) 



Here the first two terms describe the cavity mode including its 
coupling to the external continuum; the third term is the cou- 
pling between the cavity mode and the mechanical oscillator; 
Go = oo r /L c is the coupling constant with L c the cavity length; 
the rest of the terms describes the dynamics of the effective 
oscillator (left and middle boxes in Fig. [2]), with summa- 
rizing the dynamics of the cooling mode and the thermal heat 
bath, as well as their coupling with the original mechanical 
oscillator. 

In the rotating frame at the laser frequency, the linearized 
equations of motion for the perturbed part — variation around 
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the steady-state amplitude — read: 

m (x + K m .x + <a 2 x) = F B a + F q , 
a + (K r /2 + iA)d = —iGox + y/Krd in , 

where the back-action force Fqa is defined as: 

Fba = -hGo(a + cf), 



(23) 
(24) 



(25) 



and we introduce Go = aGo with a being the steady-state am- 
plitude of the cavity mode and a m is the annihilation operator 
of the input vacuum field. The cavity output a out is related to 
the cavity mode by: 



K r a. 



(26) 



with K r the decay rate (the bandwidth) of the readout mode. 

In the steady state, these equations of motion can be solved 
more easily in the frequency domain. Starting from the me- 
chanical displacement, we get 



x(co)=x(co)[F BA (co)+F,{co)}. 



(27) 



Here we have ignored modification to the mechanical re- 
sponse function % due to the readout mode — a term propor- 
tional to Gq, assuming that the pumping power is low. For the 



(28) 



cavity mode, we invert Eq. ( |24| i and obtain 

Gqx(co) + iy/K^di n (<X)) 



CO 



■A+iK r /2 



which leads to 



Fba(o)) = 



2tiG ^K r /2[(K r /2~i(o)\<i +Av 2 ] 
(0) - A + iKr/2) (CO + A + iKr/2) 



(29) 



withv! = (a in +flJ n )/\/2and v 2 = [a ia ~ a\ B ) / (^/2i) being the 
amplitude quadrature and the phase quadrature of the input 
field, which has fluctuations at the vacuum level. When com- 
bining with Eq. p6| ), we obtain the output amplitude quadra- 
ture 



Fl(co) = [a ou ,((0) +a J ont (-CO 

(A 2 - fC, 2 /4 - (0 2 )vi - K r Av 2 - 



(co - A + iKr/2) (co + A + iKr/2) 



(30) 



whose spectrum is measured experimentally. We put the 
above formula into the same format as Eq.fTJfr by normal- 
izing Y\ with respect to the mechanical displacement x, and 
introduce y(co) and the corresponding sensing noise z(co): 



y(co) = 



(A 2 - fc 2 /4 - co 2 ) v\ - K r Av 2 



y/2K r G A 

z(co)+x{co)[F BA ((o)+F q (co)}. 



+x(co) 



(31) 



Taking single-sided symmetrized spectral density of y (see 
Appendix [A]), we obtain 



where 

S zz (co) 
S z f(co) 

Here we have used 



(A 2 -k- 2 /4-0) 2 ) 2 + k- 2 A 2 

2fc r G 2 ) A 2 
h(K r /2- ico) 



2h A Gf j K r (Kf/A + co z + A l ) 
(A 2 -k- 2 /4-0) 2 ) 2 + k- 2 A 2 



(33) 
(34) 
(35) 



(0|v J (a))vJ(a)')|0> sym = n8 ]k 8(co-co') (j,k=l,2)- (36) 

Indeed, only S z j7 depends on the sign of detuning and con- 
tributes to the asymmetry. In the resolved-sideband case 
K r <^ co m and choosing detuning |A| = co m , S z f can be ap- 
proximated as the one shown in Eq. ( fT8| ). For a weak readout 
beam, we can ignore Spp which is proportional to Gq, the out- 
put spectra around co m for the positive- and negative-detuning 
cases can be approximated as 



Svv(G))|a=±»„, 



Kr | ^K m (2(n) + lTl) 

2G 2 2mco m [(co-co m ) 2 + (K m /2) 2 } 



(37) 

As we can see, the contribution to output spectra from the 
quantum correlation has the same magnitude as the zero-point 
fluctuation of the mechanical oscillator, with a sign depending 
on the detuning. One can then obtain the dependence of the 
asymmetry factor 77 on (n) shown in Eq. |9]). 

Interestingly, even if the quantum back-action term Sfp is 
much smaller than S q FF and has been ignored, given the weak 
readout mode used in the experiment, the asymmetry induced 
by quantum correlation is always visible as long as (n) is 
small. In addition, any optical loss in the readout mode only 
contributes a constant noise background — that is symmetric 
with respect to detuning — to the overall spectrum; therefore, 
the asymmetry is very robust against optical loss, and it can 
be observed without a quantum-limited readout mode which 
is the case in the experiment. 



D. Connection with the scattering picture 

In the above, we have been emphasizing the viewpoint of 
position measurement and interpreting the asymmetry as due 
to the quantum correlation between the sensing noise and the 
back-action noise. Here we would like to show the connection 
between this viewpoint and the scattering picture in Ref. E4l 
that focuses on the photon-phonon coupling, and in addition, 
show how spectra of the amplitude quadrature measured in the 
experiment are related to emission spectra of the (anti-)Stokes 
photons that would have been obtained if we instead take a 
photon-counting measurement. 

To illustrate these, we introduce the annihilation operator b 
for the phonon through the standard definition: 



S yy (co) = S zz + 2Re[x*S zF ] + \x\ 2 [S B ff + S FF ] 



(32) 



v / h/(2mco m )(b + b" r ). 



(38) 
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and it satisfies the commutator relation: [b, b'] = 1. In the ro- 
tating frame at the laser frequency, the Hamiltonian in Eq. (j22j 
after linearization is given by: 

^ = hAa^a + 3@ Kr + hg (d + a f ) (b + tf) + hco m Pb + jfi^ , 

(39) 

where go = Go\/h/ (2mco m ). The third term is the photon- 
phonon coupling: d y b describes the anti-Stokes process — 
the absorption of a phonon is accompanied by emission of a 
higher-frequency photon; a^b' describes the Stokes process 
— the emission of a phonon is accompanied by emission of 
a lower-frequency photon. The photon emission rate of these 
two processes can be estimated by using the Fermi 's golden 
rule. Specifically, taking into account the finite bandwidth for 
the photon and phonon due to coupling to the continuum, the 
emission rate of the anti-Stokes photon at C0[ r + CO reads: 



f J dTe imr @((Q)(P(T)b(0)) 
glK m {n)Si{(o) 



(co-co m ) 2 + (K m /2) 2 ' 
the emission rate of the Stokes photon at C0i r — co reads: 

9{-(0){b{x)P{<S)) 
l)^(-co) 



(40) 



r s (a>)=f JdTe->™ 



(CO-CO m ) 2 + (K m /2) 



2 ■ 



(41) 



Here the density of state for the photons is determined by the 
cavity decay rate and detuning: 



$((£>) 



Kr/2 



(co-A) 2 + (K r /2) 2 



(42) 



Were the cavity bandwidth much larger than the mechanical 
frequency co m , the density of state &(co) would become flat 
for frequencies around ±0) m , and we would effectively have 
a scenario that is similar to the free-space Raman scattering 
as in those spectroscopic measurements of crystals |3Q| . By 
making a photon-counting-type measurement of the emitted 
(anti-)Stokes photons, one could observe an asymmetric spec- 
trum with two peaks (sidebands) around co r ± CO m of which 
the profiles are given by the above emission rates. This is also 
the case for those emission and absorption spectroscopic mea- 
surements in the ions/atoms trapping experiments ||25] - E8l . 

The situation of our experiment is however different from 
the usual free-space Raman scattering spectroscopic measure- 
ment by the following two aspects: (i) we are operating in 
the resolved-sideband regime where the cavity bandwidth is 
much smaller than the mechanical frequency and the photon 
density of state is highly asymmetric for positive and negative 
sideband frequencies depending on the detuning. This basi- 
cally dictates that we cannot measure two sidebands simulta- 
neously, and we have to take two separate spectra by tuning 
the laser frequency. In the positive-detuning case A = co,„, 
the anti-Stokes sideband is enhanced while the Stokes side- 
band is highly suppressed, as the photon density of state is 
peak around co = CO m ; while in the negative-detuning case 



A = —co m , the situation for these two sidebands swaps; (ii) 
we are using heterodyne detection scheme instead of photon 
counting, where the outgoing field is mixed with a large co- 
herent optical field (reference light) before the photodetector, 
to measure the output amplitude quadrature, and the signal 
is linear proportional to the position of the oscillator, as we 
mentioned earlier. Interestingly, there is a direct connection 
between the spectra of amplitude quadrature measured in the 
experiment and the photon emission spectra that are obtained 
if making photon-counting measurements. To show this con- 
nection, we use the fact that 

[fi(fl)),fi t (fl)')]=0 (43) 

which is a direct consequence of \y{t),y(t')] = (y is equal 
to Y\ normalized with respect to the mechanical displacement 
[cf. Eq. pi)]), and we have 



= ~[{&out(-Co')al ut {-(0)) + (al ut (co')a out {co))} . 



(44) 



Take the positive-detuning case A = co m for instance, 
flout(~ ©) contains mostly vacuum and negligible sideband 
signal due to suppression of the Stokes sideband around C0\ r — 
co m by the cavity, namely, (d out (— co')dl ut (~ co)) ps 2k8(co — 
co'). The second term gives the emission spectrum for the 
output photons shown in Eq. ( |40] i; therefore, the single-sided 
spectral density of the output amplitude quadrature reads: 

S yiF » = l+2r AS (a)). (45) 

By normalizing the spectrum with respect to the mechanical 
displacement, we have 



S } y(C0)\ A=0 



2 -[i+2r AS (o))]. 

2O 



(46) 



Similarly, by following the same line of thought, we get 

S n .(co)\A=- mm = ^[l+2T s (co)}. (47) 
2G 

The above two equations give identical results to Eq. ( (37] >. 
Therefore, the output spectra obtained in our heterodyne de- 
tection differ from those in the photon-counting measurement 
only by a constant noise floor, which originates from vacuum 
fluctuation of the amplitude quadrature. After subtracting this 
noise floor, we simply recover the emission spectra obtained 
from taking photon-counting measurement. 



GENERAL LINEAR MEASUREMENTS OF THE 
ZERO-POINT FLUCTUATION 



III. 



Based on the analysis of the specific experiment of Ref. 112411 
in the previous section, here we comment on the general fea- 
tures of linear quantum measurements involving reading out 
zero-point fluctuation of a mechanical oscillator. We start 
from discussing nature of the zero-point mechanical fluctu- 
ation in Sec.|III A proceed to discussion of measurements of 
it in Sec. Ill B and finally end in Sec. Ill C which discusses 
its effect on sensitivity for measuring external forces and the 
connection to the SQL. 
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A. The nature of zero-point mechanical fluctuation 

First of all, let us take a closer look at the nature of the zero- 
point fluctuation of a realistic harmonic oscillator, which con- 
sists of a mechanical mode with eigenfrequency co m and finite 
decay rate K m . Suppose we initially decouple the oscillator 
from its environmental heat bath and turn on the coupling at 
t = 0. In the Heinserberg picture, the position and momentum 
of the oscillator at t > will be 



*q{t) =i&ee(0+y X(t-t')P q (t')dt', 
Pq(t)=pfr< x (t)+m^ x(t-t')F q (t')dt' : 



where 



■^free (0 
Pfree(Q 

m(O m 



-K m t/2 



-K„,t/2 



i(0) cos CO m t + sin (0,„i 



mco„ 



-x(0) sin co m t + ^ cos a>„,t 



mco,. 



m K„ 



^free(0 3 



(48a) 
(48b) 

(49a) 

(49b) 



equal to h. Here we note that if K m <C 0) m , when reaching the 
steady state with if ree and qf lee decayed away, we have 



Ax q ■ Ap q 



: m(0„ 



dco 
2% 



(52) 



where Sx X = \x\ 2 S'ff- Although S$x depends on the specific 
scenario, they are all constrained by a Heisenberg-like relation 
of, 



SU(o)>2hlm X (co) 1 



(53) 



which is a straightforward consequence of the commutation 
relation in Eq. ( fTTj ). The equality is achieved at the ground 
stated This enforces the same Heisenberg Uncertainty rela- 
tion: 



&x q ■ Ap q > h/2, 



(54) 



as an ideal harmonic oscillator whose quantum fluctuations 
arise "on its own", instead of having to be driven by the sur- 
rounding environment. Therefore, in the steady state, the 
zero-point fluctuation of the mechanical oscillator can be 
viewed as being imposed by the environment due to linear- 
ity of the dynamics. 



are contributions from the free evolution of the initial 
Schrodinger operators (i.e., undisturbed by the environment), 
which decay over time, and get replaced by contributions from 
the environmental heat bath [integrals on the right-hand side 
of Eqs. ( |48a| > and ( |48b| i] . Note that for any oscillator with non- 
zero decay rate, it is essential to have bath operators entering 
over time, otherwise the commutation relation between posi- 
tion and momentum: 



[x q {t),Pq{t)]=ih 
will not hold at t > because of 

[*free(0;Pfree(0] = ihe~ 



(50) 



(51) 



This dictates that the heat bath must be such that the additional 
commutator from terms containing F q exactly compensate for 
the decay in Eq. ( |5Tj ), which leads to the quantum fluctuation- 
dissipation theorem (see e.g., Ref. l37l ). 

It is interesting to note that this "replenishing" of commu- 
tators has a classical counterpart, since commutators are after 
all proportional to the classical Poisson Bracket. More specif- 
ically, for a classical oscillator with decay, we can write a sim- 
ilar relation for Poisson Brackets among the position and mo- 
mentum of the oscillator, plus environmental degrees of free- 
dom. The replenishing of the position-momentum Poisson 
Bracket by environmental ones, in classical mechanics, can 
also be viewed as a consequence of the conservation of phase- 
space volume, following the Liouville Theorem. A decaying 
oscillator's phase-space volume will shrink, and violate the 
Liouville Theorem — unless additional phase-space volume 
from the environmental degrees of freedom is introduced. 

Nevertheless, the definitive quantum feature in our situation 
is a fundamental scale in the volume of phase space, which is 



B. Measuring the zero-point fluctuation 

Having clarified the nature of quantum zero-point fluctua- 
tions of a mechanical oscillator in the steady state, let us argue 
that the effects seen in Ref. [24] are actually generic when one 
tries to probe such fluctuations, namely: the correlation be- 
tween sensing and back-action noise can be at the level of the 
zero-point fluctuation itself. 

Let us start our discussion here from Eq. |5]), namely, 



(55) 



and the fact that y consists of sensing noise, back-action noise, 
and finally the zero-point fluctuation of the mechanical oscil- 
lator [cf. Eq. (|8]l]: 



m = 



m 

a 



« / Z(f-T)F BA (T)dT-|-i 9 (f). 



(56) 



Here we have added a factor a, which labels the scaling of 
each term as the measurement strength which is proportional 
to the square root of the readout beam power. Let us assume 
that Eq. ( |55] l continues to hold for the same set of z and Fba, 
for a large set of a and basically the measuring device 
works for different mechanical oscillators with different mea- 
suring strength. 

Since Eq. ( |55] > remains valid for all values of a, we extract 
terms with different powers of scaling, and obtain 



[z{t), z{t')] = [Fba(0^ba(O] =0, 



(57) 



• A generalization of this to thermal states will be the fluctuation-dissipation 
theorem |37|. 
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and 



X(t'-x) [z(f), J F BA (T)]dT 

Z(f-T)[£(f'),^BA(T)]dT 
Xq(t),Xg(t')]=0, Vf,f'. 



(58) 



This becomes 

r+oo 







X(r) [C zF (t - t) C z f ( t - r)]dT - -ihx(t) , (59) 



for all values of f , where we have defined 

C zF (t'-t)=[z(t),P BA (t')] 



(60) 



Here the dependence is only through t' — t because the system 
assumed to be time-invariant. We also note that since z is 
an out-going field, C z p(t' — t) must vanish when t' — t > 0, 
otherwise any generalized force applied on the out-going field 
z(t) — detached from the mechanical oscillator — can still 
dynamically influence the mechanical motion at later times 
(future) through fBA(f'), which violates the causality lfTTl[35l . 
As proved in the App.[B] in order for Eq. ( f59] l to be satisfied 
for all possible response functions of the oscillator, we must 
have 



C zF (t) = -ih8-{t), 



(61) 



where 5_(f) is the Dirac delta function with support only for 
t < 0. In other words, 



[z(t),P BA (t')] = -ihS-(t'-t) 



(62) 



Eq. ( |56") , plus the commutation relations in Eqs. ( |57) and 
([52}, then provides a general description of linear measur- 
ing devices which do not modify the dynamics of the me- 
chanical oscillator — simply from the requirement that the 
out-going field operators at different times must commute 
[cf. Eq. <|55)]. In particular, non-vanishing commutator 
[x q (t), x q (t')], which underlies the existence of the zero-point 
fluctuation, is canceled in a simple way by the non-vanishing 
commutator between the sensing noise and the back-action 
noise [cf. Eq. ([62])], 

Now turn to the noise content of the output y(t), i.e., the 
spectrum, 



Syy 



^+2Re[x*S Z F} + a 2 S FF 



(63) 



Let us consider experiments with relatively low measurement 
strength, so that the first term S zz /a 2 from the sensing noise 
dominates the output noise. The next-order terms contain: 
(i) correlation between the sensing noise and the back-action 
noise — S z f', and (ii) the mechanical fluctuation — Sf x - If we 
assume nearly ground state for the mechanical oscillator 



which, for K m <C CO„ 



gives 
da> 

2jZ 



2m(O m 



(64) 



(65) 



If S z f(co) does not change noticeably within the mechanical 
bandwidth, then 



^2Re[x* (co)S zF (co)} « -— !— ImS^c^) . (66) 
2% 2m(Q m 



Because of Eq. ( [62) , the typical magnitude for S z f is natu- 
rally 11 



\S zF 



(67) 



Therefore, contributions to the output noise from quantum 
correlation S z f and mechanical fluctuation S q xx can generically 
become comparable to each other when the mechanical os- 
cillator is approaching the quantum ground state. The result 
presented in Ref. [24| therefore illustrates two typical cases of 
this generic behavior [cf. Eq. ([20)]. 



C. Measuring external classical forces in presence of 
zero-point fluctuation 

Finally, let us discuss the role of zero-point fluctuation in 
force measurement, when the mechanical oscillator is used 
as a probe of external classical forces not far away from the 
mechanical resonant frequency. The force sensitivity of such a 
linear measurement device, in terms of spectral density Sf, is 
obtained by normalizing the displacement sensitivity S yy with 
respect to the mechanical response function %: Sf = S yy /\x\ 2 . 
Specifically, from Eq. ( fl3) , we have 



\xi.o))\ 2 



2Re 



S z f((Q) 
X{®) 



+ S B FF ( ( o) + S FF ((o) 



(68) 



Because of the commutation relations in Eqs. |57) and ( |62) , 
a Heisenberg Uncertainty Relation exists among the spectral 
densities of z and Fra, and that is 



S zz ((o)S FF {co) - S zF {co)S Fz (co) > h z 



(69) 



When the the sensing noise z and the back-action noise Fba 
are not correlated — S z f = Sf z = 0, we have 



S zz ((o)S FF (co)>h z 



(70) 



The above inequality represents a trade-off between sensing 
noise z and back-action noise Fba- Correspondingly, the force 
sensitivity will have a lower bound : 

S F ((D)\ SzF=0 =-^^+S P F((D)+^ F (a ) ) 
\XW\ 
2h 

> rT-^y + (4<n> + 2)hmK m (O m . (71) 



' In general, the commutator does not impose any bound on the cross corre- 
lation. Here, in a strict sense, is an order-of-magnitude estimate. 
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FIG. 3: (color online) Figure illustrating that total quantum limita- 
tion 5^ tot (red) for force sensitivity and contribution from zero-point 
fluctuation S Z p (blue). For clarity, we divide both by the SQL and use 
log-log scale. 



c Qtot2 




10re m Z^y^^ 


s z / 2 






\ /cQtotl 
\ / 









FIG. 4: (color online) Figure illustrating effect of the mechanical 
decay rate (bandwidth) K m on 5^ tot (solid) and S F P (dashed) — the 
larger the mechanical bandwidth, the lower the force sensitivity (This 
plot is also in the log-log scale). 



If the mechanical oscillator is in its quantum ground state, 
namely (n) = 0, we obtain: 

S F {C0)> + 2hmK m (O m = Sf 0t . (72) 

The first term is the usual Standard Quantum Limit (SQL) for 
force sensitivity with mechanical probes ifTTl fT2l : 

S S F QL = = 2^(0)2-^)2 + ^0)2. (73) 

The second term, 

Sf = 2hmK m (0 m , (74) 

arising from the zero-point fluctuation due to mechanical 
quantization, also limits the sensitivity. As we can learn from 
Eqs. ((68), ( |69) and ( |72) , the quantum limit, can be surpassed, 
in principle indefinitely, by building up quantum correlations 
between the sensing noise z and the back-action noise Fba — 
in practice the beating factor will be limited by the available 
optical power and the level of optical losses. However, the 
limit imposed by zero-point fluctuation cannot be surpassed, 
and it can only be mitigated by lowering K m , i.e., increasing 
the mechanical quality factor. 

Braginsky et al. Il38l argued that mechanical quantization 
does not influence the force sensitivity when measuring clas- 
sical forces with mechanical probes — one only needs to eval- 
uate the quantum noise due to the readout field. But these 
authors had specifically pointed out that they were focusing 
on ideal mechanical probes with infinitely narrow bandwidth 
(K m — > 0) and observations outside of that frequency band. 
This is close to the actual situation of free-mass gravitational- 
wave detectors, in which the mechanical oscillator is the dif- 
ferential mode of four mirror-endowed test masses hung as 
pendulum with eigenfrequencies around 1 Hz and very high 
quality factor, while the detection band is above 10 Hz, well 
outside the mechanical resonance. Indeed, from Eqs. ( f73) and 
( |74) , we see that the effect of zero-point fluctuation is only 
significant not far away from resonance — which confirms 
Braginsky et a/.'s result. More specifically, if K m <C C0 m , we 



can write, for | CO — C0 m | <C CO,,,, 

In particular, the limit imposed by zero-point fluctuation is 
equal to SQL on resonance, and becomes less important as 
| co — co„, | becomes comparable to or larger than the half band- 
width rc,„/2, as illustrated in Fig. [5] Note that on an absolute 
scale: S F ® L (co) is lower near the mechanical resonance, while 
Sp((o) is independent from frequency; at any frequency, low- 
ering Km, while fixing CO,,, and keeping the oscillator at ground 
state, always results in lower noise, as illustrated in Fig. [4] 
Suppose we are free to choose from ground-state mechanical 
oscillators with different co,„ and K m as our probe, and that we 
are always able to reach the SQL at all frequencies, then: (i) 
if we know the frequency content of target signals, we can 
choose probes that are closely resonant with the target, and 
(ii) regardless of signal frequency, probes with lower K m , or 
equivalently, higher mechanical quality factor, always provide 
better force sensitivity. 

IV. CONCLUSION 

We have shown, within the framework of quantum mea- 
surement theory, that the asymmetry in output spectra ob- 
served by Safavi-Naeini et al. can be explained as due to the 
quantum correlation between the sensing noise and the quan- 
tum back-action noise; this experiment therefore provides a 
clear signature of quantum back-action onto mechanical sys- 
tems. More broadly, we have shown that having quantum- 
noise correlations showing up at the same level as the zero- 
point fluctuations is a generic feature of measurements that 
attempt to measure the zero-point fluctuation. We have further 
shown that when an experimentally prepared ground-state me- 
chanical oscillator is used as a probe for classical forces near 
its resonant frequency, its mechanical quantization — through 
zero-point displacement fluctuation — does impose an addi- 
tion noise background. This additional noise vanishes only 
if the oscillator's bandwidth approaches zero, i.e., when the 
oscillator becomes ideal. 
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Appendix A: Symmetrized cross spectral density 

In this article, as in Ref. Ifl4ll . we use the single-sided sym- 
metrized cross spectral density, which, given a quantum state 
| y/), is defined between a pair of operators A and B as: 



S AB {co)8{(o - a') = -(xi/\A{(o)B\(o')\\i/) sym 



In 



(\l/\A(co)B T {a')+B T ((o')A((D)\\if). (Al) 



The symmetrization process here allows us to preserves bilin- 
earity of S on its entries, i.e., 



Sa,c 1 B+c 2 C — c\Sab + C 2^AC : 

S Ci a+c 2 b,c = c i Sac + ciSbc ■ 
More importantly, we can show that 

Saa >0 



(A2a) 
(A2b) 



(A3) 



for any field A, even if [A(o), A + (co')] ^ 0. The positiv- 
ity ( |A3[ > allows us to interpret Saa as the fluctuation variance 
per unit frequency band — as in the classical case. 



Appendix B: Commutation relation between z and F 

Defining 

f{t)=C zF {t) + ih8-{t) (Bl) 
we convert Eq. |59) into 

°%{x)\f{t-x)-f{-t-x)]dx = 0. (B2) 



Assuming analyticity of the Fourier Transform of /(f), it must 
be written as 



y Ih 



(B3) 



with a>k all located on the upper half of the complex plane (not 
including the real axis). Fourier transforming Eq. \B2\ gives 
us 



k 



h ft 

co — Ok a>-Q>t 



= 0, toe 



(B4) 



Because the set {cty} is within the upper-half complex plane 
(excluding the real axis), the set {co^} must be within the 
lower-half complex plane (excluding the real axis) — and the 
two sets do not intersect. For this reason, Eq. ( |B4| i requires /. 
to all vanish, and hence 



C zF (t) = -ih8-(t). 



(B5) 
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